home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / EDUCNOMY / ASTROSET.LZH / STAY.BAS < prev    next >
BASIC Source File  |  1985-12-20  |  896b  |  42 lines
  1. 10 REM   LENGTH OF STAY
  2. 12 REM   INSIDE 1 A.U.
  3. 14 REM
  4. 16 K=0.01720209895
  5. 18 P=3.1415926536
  6. 20 INPUT "PERIHELION DIST. ";Q
  7. 22 INPUT "ECCENTRICITY     ";E
  8. 24 PRINT
  9. 26 T=0
  10. 28 IF E>1 OR E<0 THEN 22
  11. 30 IF Q>=1 THEN 78
  12. 32 IF Q<0.0046524 THEN 82
  13. 34 IF E<1 THEN 54
  14. 36 REM
  15. 38 REM   PARABOLIC CASE
  16. 40 REM
  17. 42 T2=1/Q-1
  18. 44 T1=SQR(T2)
  19. 46 Q3=Q*Q*Q
  20. 48 T=SQR(2*Q3)*(3*T1+T1*T2)/(3*K)
  21. 50 GOTO 78
  22. 52 REM
  23. 54 REM  ELLIPTICAL CASE
  24. 56 REM
  25. 58 IF E=0 THEN 86
  26. 60 IF Q*(1+E)/(1-E)<=1 THEN 86
  27. 62 A=Q/(1-E)
  28. 64 C=(A-1)/(A*E)
  29. 66 IF C<>0 THEN 70
  30. 68 E1=P/2:  GOTO 76
  31. 70 S=SQR(1-C*C)
  32. 72 E1=ATN(S/C)
  33. 74 IF C<0 THEN E1=E1+P
  34. 76 T=A*SQR(A)*(E1-E*SIN(E1))/K
  35. 78 PRINT "DAYS: ";2*T
  36. 80 GOTO 88
  37. 82 PRINT "COLLISION WITH SUN"
  38. 84 GOTO 88
  39. 86 PRINT "ALWAYS INSIDE 1 A.U."
  40. 88 END
  41.  
  42.